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1 Introduction 

Recent data on hadron production in central nucleus-nucleus (A+A) collisions at the 
CERN SPS are compatible with the hypothesis that a Quark Gluon Plasma J]] is created 
in the early stage of the interaction 0. 

This interpretation however requires the assumption that the produced matter is close 
to thermal and chemical equilibrium. Thus it is important to measure the level of equili- 
bration reached in nuclear collisions. 

The method to study event-by-event fluctuations of kinematical variables ("thermal" 
fluctuations, TF-method) was proposed in 1992 || and recently used for the analysis of 
central Pb+Pb collisions at 158 A-GeV by the NA49 Collaboration ||. First results indi- 
cate that the analysis of event-by-event fluctuations adds crucial information concerning 
the dynamics of A+A collisions and in particular allows to reject purely non-equilibrium 
approches, like initial state scattering models ||. The amount of "thermal" fluctuations 
was recently calculated for matter in equilibrium ||. 

In this paper we propose a method to study event-by-event fluctuations of the particle 
composition of the final state of high energy collisions ("chemical" fluctuations, CF- 
method). We expect that this method will allow to determine whether or not chemical 
equilibrium is reached in high energy collisions. As it is based on the TF-method, we start 
(Section 1) from a brief description of the basic idea and the formalism of the TF-method 
@. In Section 3 we introduce the CF-method. A simple numerical example is presented 
in Section 4. 



2 Description of the TF-Method 

Let us suppose that A+A collisions can be modeled as a sum of independent nucleon- 
nucleon (N+N) interactions. In this case event-by-event fluctuations in A+A collisions 
are given by a superposition of the fluctuations present in N+N interactions. Additional 
fluctuations are introduced when the number of superimposed N+N interactions varies 
from event to event due to e.g. variation in the collision geometry. 

A statistical method which allows to remove the influence of trivial geometrical fluc- 
tuations and the effect of averaging over many particle sources was introduced in Ref. ||. 
It was proposed to quantify the fluctuations by the so-called $ variable. As an example 
we consider here a construction of $ for the case of transverse momentum fluctuations. 

For every particle i one defines: 

Zi = Pt x ~ Pt, 

where pr is the mean transverse momentum calculated for all particles from all events 
(the inclusive mean). Using Zj one calculates for every event 

N 
i=l 



where N is the number of analyzed particles in the event. The fluctuation measure^], <3> PT , 
is then defined as: 
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where (N) and (Z 2 ) are averages (of event-by-event observables) over all events and the 
second term in the r.h.s. is the square root of the second moment of the inclusive z 
distribution. 

By construction ||, the <3> PT value for A+A collisions is equal to the $ PT value for N+N 
interactions in the case in which A+A collisions are pictured as a sum of independent N+N 
interactions. If the particles are produced independently the value of <3> PT is equal to zero. 

3 Formulation of the CF-Method 

The TF-method can be converted into a method allowing to study event-by-event fluc- 
tuations of the "chemical" composition (relative number of different hadronic states) of 
the collisions (CF-method). In the latter method the basic quantity in which fluctuations 
are analyzed is the number of particles of a given type produced in a single collision. This 
number substitutes the transverse momentum of particles used in the TF-method (see 
Section 2). 

A formal trick which allows for a direct conversion of the TF-method to the CF- 
method is based on the substitution of the kinematical variable used for "thermal" fluc- 
tuation studies (e.g. pr as discussed in Sect. 2) by a function defined as: 

'(*'■»•) = {S ; iv* (2) 

where i is a particle index, h l is the particle type of the particle % and ho is the particle 
type selected for the fluctuation analysis. Then the variable z takes the form: 



Zi = S(h\h ) -5(h,ho), 



where the 8(h, ho) is the mean value of 5(h l , h ) calculated for all particles from all events 
(the inclusive mean). 5(h, h ) gives the probability that any particle of the event ensamble 
is of the type ho. For every event the sum: 

N 
i=l 

is calculated where the summation runs over all iV particles. Finally the "chemical" 
fluctuation measure, Q(h ), can be calculated as: 
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2 In the original paper || the $ PT is called AD. Here we follow notation introduced by the NA49 
Collaboration, which relates the name of the fluctuation measure to the variable in which fluctuations 
are studied. 



where (N) and (Z 2 ) are averages (of event-by-event observables) over all events and the 
second term in the r.h.s. is the square root of the second moment of the inclusive z 
distribution. 

As follows from the construction the important statistical features of the $(/io) variable 
are identical to the features of the $ variable used in the TF-method (e.g. $ PT ). We list 
them below. 

• For a system which is an independent sum of elementary processes the value of $(/^o) 
is equal to the value of &(h ) calculated for a single elementary process independent 
of the number of superimposed elementary processes and its distribution in the 
analyzed event sample. 

• In the model in which particles are produced independently from each other the 
value of Q(h ) is equal to zero. 

4 Numerical Example 

We first note that particle production in elementary processes (e.g. p+p interactions) is 
correlated not only in momentum space |7j] but also when the "chemical" composition is 
considered. In order to illustrate this statement by experimental data || ||| we show in 
Fig. 1 the ratio of the mean K s multiplicity to the multiplicity of negatively charged 
hadrons, n~, as a function of n~ for p+p interactions at 200 GeV/c. The data show 
that the K s multiplicity decreases significantly with the event multiplicity. A similar 
correlation can be expected for K~ multiplicity, and in fact, it is observed in string models 
of p+p interactions jlQf . It means that the probability that a negatively charged hadron 
is K~ meson depends on the multiplicity of negatively charged hadrons in the event. 
Thus particles are not produced independently and therefore we expect that the value 
of $(fr~) calculated for negatively charged hadrons for p+p interactions is not equal to 
zero. In order to make numerical estimation of the effect we use a simple parametrization 
of hadron production in p+p interactions at 200 GeV/c. 

It is assumed that a dependence of (K~)/n~ on n_ is similar to the dependence of 
(K s )/n_ on n_ shown in Fig. 1. The multiplicity distribution of negatively charged 
hadrons is calculated using the parametrization from Ref. |lTJ. Further we assume that 



the multiplicity distribution of K mesons for a fixed multiplicity n_ is given by the 
binominal distribution i.e.: 



%«-)= -\P K {n-) nK {l-PK{n-)) n - nK , (4) 




where nx is the kaon multiplicity and Pxin-) is the probability that a negatively charged 
hadron is if "-meson. Following the data presented in Fig. 1 the Pft-(n_) is parametrized 
as: 

P ,„ n_J Pa:(1) - (n- - 1) ■ A : P K (1) > (n_ - 1) • A 

P ^ n -)-{ o : P K (1) < (n_ - 1) • A (5) 

In order to study the dependence of <&(K~) on the correlation between P^nS) and n_ 
the parameter A is varied between (no correlation) and 0.02 (a correlation as suggested 



by the data presented in Fig. 1). The parametrizations of Px{n-) obtained for these two 
extreme values of A are indicated by dashed lines in Fig. 1 for Pr-(I) = 0.13. 

The dependence of <&(K~) on the parameter A is shown in Fig. 2 for two values of the 
parameter Pr-(I): 0.13 (solid line) and 0.26 (dashed line). For A = the value of <&(K~) 
is equal to zero, as expected from the definition of $. The value of Q(K~) increases 
with increasing A, i.e. with the increasing correlation between K~ yield and the event 
multiplicity. 

Using the above model (P^(l) = 0.13 and A = 0.02) we estimated that the number 
of events needed to obtain 10% statistical error on $(iT~) is about 10 5 . 

5 Summary 

We proposed a method to study event-by-event fluctuations of the "chemical" composi- 
tion of the final state of high energy collisions. A simple numerical example of fluctuations 
of the number of i^~-mesons in the sample of negatively charged hadrons was considerd. 
The method can be used to study a change of the magnitude of the "chemical" fluctua- 
tions when changing the size of the colliding systems (p+p, p+A and A+A) and/or when 
changing the collision energy. 

We expect that the method will allow to determine whether or not chemical equili- 
bration is reached in high energy nucleus-nucleus collisions. 
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Figure 1: The dependence of the ratio {Ks)/n_ on negative hadron multiplicity, n_, for 
p+p interactions at ,/s = 20 GeV. 
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Figure 2: The dependence of the fluctuation measure Q(K ) on the paremeter A (see 
text) for two values of the parameter Px-(l): 0.13 (solid line) and 0.26 (dashed line). 



